Equilibrium topology of the intermediate state in type-I superconductors of different 

shapes 
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High-resolution magneto-optical technique was used to analyze flux patterns in the intermediate 
state of bulk Pb samples of various shapes - cones, hemispheres and discs. Combined with the 
measurements of macroscopic magnetization these results allowed studying the effect of bulk pinning 
and geometric barrier on the equilibrium structure of the intermediate state. Zero-bulk pinning discs 
and slabs show hysteretic behavior due to topological hysteresis - flux tubes on penetration and 
lamellae on flux exit. (Hemi)spheres and cones do not have geometric barrier and show no hysteresis 
with flux tubes dominating the intermediate field region in both regimes. It is concluded that flux 
tubes represent the equilibrium topology of the intermediate state and that the laminar structure 
is unstable towards Lorentz or condensation energy forces. Real-time video is available in [24[. 



Pattern formation in strongly correlated systems is a 
topic of incessant interest for a broad scientific commu- 
nity At a first sight, type-I superconductors repre- 
sent a perfect physical system where it is relatively easy 
to tune the parameters and try to understand the physics 
behind observed topology of the intermediate state. In 
fact, analogies to type-I superconductors extend from as- 
trophysics 2] to the physics of ice 0] . The fundamental 
problem, however, is that in a finite system it is impossi- 
ble to predict 2D and moreover 3D pattern based solely 
on the energy minimization arguments 0]. The pattern 
has to be assumed and then its geometrical parameters 
are determined from the minimization. Back in the 1930s 
Lev Landau suggested a simple stripe model, which was 
possible to analyze analytically [f], [y] . Later refinements 
(such as domain widening and/or branching) tried to ad- 
dress apparent inconsistencies between the model and the 
experiment @> B @|- Still a comprehensive description 
has never been achieved with the main problem being 
multiple observations of the closed-topology structures 
(flux tubes) in best samples. Experimental and theoret- 
ical effort is growing to obtain general understanding of 
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the problem 0, El El • 

Here we outline several fac- 
tors that influence and often determine the topology of 
the intermediate state and must be taken into account 
by a successful theory. 

First issue is (bulk) flux pinning. In the overwhelm- 
ing majority of published papers this question is simply 
omitted. We have shown that tubular topology is fairly 
robust, but is destroyed and transformed into a laminar 
pattern by the structural defects. In more disordered 
samples with significant bulk pinning, a non-equilibrium 
dendrite - like topology of the intermediate state is ob- 
served [lij • Recently, it was shown that close to H c lam- 
inar structure transforms into tubular upon applying a 
small-amplitude AC field to shake it out of the metastable 
state [ltj. Fortunately, it is easy to distinguish between 
pinning-induced and topological hysteresis. The former 
becomes larger at the lower fields and reaches maximum 
at H = 0. The latter is maximal around H = H c (1 — N) 
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FIG. 1: M(H) loop measured in a Pb single crystal at T = 4.5 
K. Larger open symbols show measurements each after field 
cooling at indicated field to 4.5 K. Small open symbold show 
minor hysteresis loops. 



(where N is the demagnetization factor) and should van- 
ish in the limit of H = 0. We note that a concept of 
topological hysteresis as applied to type-I superconduc- 
tors was introduced to describe the irreversibility in a 
macroscopic response due to different topologies of the 
intermediate state 



e.g., sample magnetic moment 14 1. 
It was later used to describe a subtle specific issue of a 
crossover between tubes and laminae in films near the H c 

Second crucial ingredient to understand the interme- 
diate state is sample shape and geometry. Films are rel- 
atively easy to make and work with, but they provide 
only very limited experimental information and are at 
the borderline of applicability of most theories. The ef- 
fective penetration depth depends on film thickness and, 
in case of Pb, the sample exhibits behavior of a type- 
II superconductor below about 0.1 fim The phase 
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interior The barrier is different for flux exit, there- 
fore tubes were considered to be a result of the metastable 
state determined by the geometric barrier. 

In this paper we study pinning-free samples of vari- 
ous shapes - discs, (hemi) spheres and cones. The latter 
two geometries do not have the geometric barrier (as di- 
rectly proven by the observations), yet show that flux 
tubes exist both upon flux entry and exit proving that 
flux tubes represent the equilibrium topology of type-I 
superconductors. 

Quantum Design MPMS magnetometer was used for 
magnetization measurements. Magneto-optical (MO) 
imaging was performed in a pumped flow-type optical 
4 He cryostat using Faraday rotation of a polarized light 
in Bi - doped iron-garnet films with in-plane magnetiza- 
tion 14l| . In all images bright regions correspond to the 



FIG. 2: Structure of the intermediate state in a disc - shaped 
Pb single crystal at 5 K. Left column - increasing magnetic 
field after ZFC. Right column - decreasing field. 



transition in the vicinity of H c is modified compared to 
the bulk case [16| and huge demagnetization factor causes 
any small defects on the film edge to act as centers of pre- 
mature flux penetration. In addition, geometric barrier 
[l^ | plays dominant role and actually determines the flux 
pattern 17[ . On the other hand, experiments with thick 
pinning-free type-I superconducting slabs and disks have 
consistently showed the tubular pattern upon flux pene- 
tration and laminar structure pattern upon flux exit. It 
was directly observed in single crystals of Sn back in 1958 
[13], Hg and Pb [IQ. The tubular pattern was re- 
produced numerically Il9| and metastability of slabs was 
theoretically analyzed [21| . Ginzburg - Landau equations 
have stable multiquanta solution for arbitrary tube size 
9]. Still, the unsatisfactory fact is that in any shape 
where there are two parallel surfaces perpendicular to 
the magnetic field, there will be a geometric barrier that 
promotes an edge instability and drives the tubes into the 



normal state and dark regions to the superconducting 
state. Due to very large volume of images and video, we 
could only include a small subset of data in this paper. 
For a complete coverage, including real-time video, see 
Ref.jii. 

We begin with a single crystal of lead in form of 
a disk of diameter d = 5 mm and thickness t = 1 
mm. Four crystals of different orientations, (110) and 
(100), and from different companies, MaTecK GmbH 
and Metal Crystals and Oxides Ltd., were studied. The 
MaTecK crystals showed lowest residual magnetic hys- 
teresis and correspondingly clearer patterns of the inter- 
mediate state. Figure [T] shows magnetization loop mea- 
sured in a (100) - oriented Pb single crystal at T = 4.5 
K. The hysteresis vanishes at H — > and minor hystere- 
sis loops (shown by smaller open symbols with the field 
sweep direction indicated by arrows) show no hystere- 
sis. Larger open symbols show result of a field-cooling 
experiment. They coincide with the data obtained by 
sweeping magnetic field down. These observations im- 
plies zero bulk pinning and we assert that the hysteresis 
comes from the difference in topologies of the intermedi- 
ate state upon flux entry and exit. 

This assertion is directly confirmed by the magneto- 
optical images shown in Fig. [2j Left column shows 
flux penetration, right column - flux exit. There is ob- 
vious difference between the topologies of the flux pat- 
terns. The tubes have sizes from [im at low fields to 
sub-mm at higher fields. We note that tubular structure 
is only observable in samples with no pinning (pinning 
leads to dendrites [3). Online vide-figures also show 
robustness of the flux tubes upon penetration of a tilted 
field [HI • It is important to note that although lami- 
nar structure appears on a hysteretic branch, it matches 
the field-cooled data indicating that this irreversibility 
is not due to macroscopic flux gradient and cannot be 
"quenched" by the annealing. We also note that tubu- 
lar pattern appears on an ascending branch that behaves 
as a textbook M (H) loop, providing additional evidence 
for its equilibrium state. 
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FIG. 3: Magnetisation loop in a "perfect superconducting" 
sphere (solid symbols) and a hemisphere (open symbols) at 
T = 4.5K. 




In fact, observation of a closed topology on flux en- 
try and open on flux exit is quite typical for any sample 
with two flat surfaces perpendicular to the applied field. 
The flux structure is governed by the geometric barrier 
[ill]. When a magnetic field is increased, Meissner cur- 
rents flow on both surfaces even in a thick disc 25]. Any 
closed-flux object will be instantaneously driven to the 
sample center by the Lorentz force. Therefore, it will 
appear as if the flux tubes pile up from the center out- 
wards, exactly as it is observed in Fig. [2] To eliminate 
the geometric barrier the sample should be in form of an 
ellipsoid. In that case, the Lorentz force is exactly bal- 
anced by the condensation energy force upon flux pen- 
etration. A comprehensive study of the shape effect on 
magnetic hysteresis was reported in Ref. ]22j where the 
authors arrived to a clear conclusion that shape plays an 
important role and showed that ellipsoidal sample has no 
hysteresis. Figure [3] shows experimentally such "perfect" 
M (H) loops measured in a Pb sphere (solid symbols) and 
a hemisphere (open symbols). The sphere was produced 
by dropping molten lead (99.9999% pure) in an inert at- 
mosphere. The hemisphere was cast into a copper mould 
in an inert atmosphere and subsequently polished and 
annealed in vacuum at 250 °C for 24 hours. Still, the 
surface was not perfect and defects are seen in the imag- 
ing. Figure [4] shows the results for a hemisphere. (Inset 
shows another hemisphere sample). The major result is 
evident - the geometric barrier is no longer present (no 
dome-like formation of flux in the center) and flux tubes 
appear both ways - on flux entry and flux exit. Real - 



FIG. 4: Flux penetration (left column) and exit (right col- 
umn) in a Pb hemisphere (d = 4 mm) at T — 4 K and 
indicated magnetic fields. Inset shows tubular pattern in a 
different hemisphere sample. 



time video of another hemisphere is available at [24|. 

The video also shows flux tubes formation upon field 
cooling as well as warming up from the pure Meissner 
state. Therefore, in a hemisphere flux tubes appear from 
all four possible ways to reach a particular point inside 
the intermediate state domain on an H — T phase dia- 
gram. 

Another shape where geometric barrier is not present 
is a cone. The cones are interesting, because one can go 
from an obtuse to an acute cone and quite possibly the 
topology of the intermediate state will change. However, 
it is very difficult to produce stress-free samples. We 
report data on the obtuse cone (4 mm diameter, 1mm 
height) in Fig. The rounded shape of the curve is 
apparently due to large local demagnetization at the cor- 
ners. Overall, the curve is quite reversible. 

Figure [5] shows magneto-optical imaging obtained in 
a conical sample. Despite presence of some defects and 
artifacts (from grease and surface imperfections), there 
is a clear absence of the geometric barrier and presence 
of flux tubes. 

We conclude that tubular structure appears to be the 
equilibrium topology of the intermediate state, because 
it is always observed on flux penetration in samples of 
any shape, as well as on flux exit and upon warming and 



4 




H (Oe) 



state at high fields, which breaks into tubes at smaller 
fields (as was foreseen by Landau @). In samples with- 
out geometric barrier, there is always a condensation en- 
ergy gradient that destroys the laminar pattern. Our 
more recent experiments also show transformation of the 
laminar pattern into tubular by applied external current. 
Importantly, our observations emphasize that different 
topologies result in actual macroscopic (topological) hys- 
teresis in magnetization as evident from FigfT] 
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53285, the Alfred P. Sloan Foundation and by the Direc- 
tor for Energy Research, Office of Basic Energy Sciences. 



FIG. 5: Magnetization loop in a cone at T — 4.5 K. 




FIG. 6: Flux penetration (left column) and exit (right col- 
umn) into a cone-shaped sample at 4 K and at indicated fields. 
Inset shows a zoom of the tubular pattern. 



cooling in constant field in samples without geometric 
barrier. Laminar structure on the other hand is unstable 
in the presence of force - either Lorentz or condensation 
energy. In flat samples, geometric barrier is present on 
penetration, but not on exit. Combined with internal 
magnetic pressure for flux to exit the latter preserves 
the laminar structure as metastable, but flux-percolative 
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